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Abstract 

We show that it is possible to construct a lattice Schrodinger func- 
tional for standard Wilson fermions, where the expectation values of 
7^-5-even operators are O(a) improved, up to terms coming from the 
boundaries. 
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1 Introduction 



It has been shown in ref. that it is possible to improve the approach to the 
continuum limit of correlation functions in lattice QCD with standard Wilson 
fermions by taking arithmetic averages (Wilson averages - WA's) of vacuum 
expectation values (v.e.v.'s) computed in theories regularized with opposite 
values of the Wilson parameter, r. Improved energies and matrix elements 
can be obtained by similarly averaging the corresponding physical quantities 
separately computed within the two regularizations. The same result can be 
obtained by replacing the WA with a linear combination of v.e.s.'s computed 
with the same value of r, but opposite values of the quark mass, m q (mass 
average - MA). In our notations m q = M — M cr , with M the bare quark 
mass parameter and M cr the critical mass. The relevant coefficient in the 
linear combination of the MA is the 71$ parity (see eq. (jH.lj) below) of the 
operator, O, whose v.e.v. one is computing. Notice that in taking the limit 
m q — > the two terms of the MA may not be equal if spontaneous chiral 
symmetry breaking occurs. 

To deal with the problems related to the spectrum of the Wilson-Dirac 
operator j2], twisted-mass lattice QCD (tm-LQCD) |3j should be better used 
for actual computations. The choice uj = ±7r/2 for the twisting angle (maxi- 
mal twist) is particularly interesting, as 0(a) improved estimates of all inter- 
esting physical quantities can be obtained even without averaging data from 
lattice formulations with opposite Wilson terms. 

These results are valid both in the quenched approximation and in the 
full theory. Indeed it was also shown in ref. |lj how one has to deal with the 
case of mass non-degenerate quark pairs to get 0(a) improvement and at the 
same time a real and positive quark determinant. 

Maximally twisted tm-LQCD can also be shown to be a sufficiently flexi- 
ble regularization scheme to allow neat solutions [13 Ej of the difficult problem 
that goes under the name of "wrong chirality mixing" which has up to 
now prevented a (reliable) lattice evaluation of the AI = 1/2 non-leptonic 
kaon decay amplitudes, if Wilson fermions are employed. By coupling the 
strategy proposed in ref. jlj for the sea quark regularization with suitable 
"twistings" of the Wilson terms of the valence quarks of various flavours, 
one can rigorously prove that it is possible to get rid of all "wrong chiral- 
ity mixings" in the evaluation of the matrix elements of the CP-conserving 
AS* = 1,2 effective weak Hamiltonian. 

Striking confirmation of the viability of the approach outlined above and 
of the remarkable properties of tm-LQCD at |u;| = tt/2 has come from the re- 
cent works of refs. 1313110], where (quenched) studies of the scaling behaviour 
of the theory were carried out down to rather small values of lattice spacing 
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and pion mass. Indeed, lattice results for pion masses as low as ~ 250 MeV 
show surprisingly small cutoff effects from a ~ .12 fm to a ~ .05 fm, if the 
optimal value of the critical mass [TTJ [12J [TBI El EH] is employed. A careful 
study of the non-trivial phase structure of the unquenched [E3 EH HZ| theory 
was also carried out in refs. fHl EH] with various choices of the gauge action 
(standard plaquette, DWB2 [20] and tree- level Symanzik improved). 

A complementary role has been played by the Schrodinger functional 
formulation |2T|l22] which has proved to be an invaluable tool in applications, 
especially because it provides a workable scheme for the computation of the 
running of the gauge coupling constant and the quark masses (2JJ ES], as 
well as for the non-perturbative evaluation of renormalization constants [21]. 
Thus a natural question to ask is whether the sort of 0(a) improvement 
one could get in the infinite volume formulation of lattice QCD can also be 
obtained in the Schrodinger functional framework. 

In this paper we construct a modified Schrodinger functional describ- 
ing the gauge interactions of a flavour doublet of massless standard Wil- 
son quarks, in which fermions are endowed with a sort of twisted boundary 
conditions. We shall show that, under the assumption that IR effects are 
screened by the natural cutoff provided by the finite extension of the time 
coordinate and up to 0(a) (fermionic and gluonic) contributions coming from 
the boundaries, the arithmetic average of the expectation values of 7?.5-even 
(multi-local, gauge invariant, multiplicative renormalizable - m.r.) opera- 
tors, O, computed with opposite values of the Wilson parameter (WA's) are 
O(a) improved. Like in the infinite volume case, the symmetries of the theory 
allow to conclude that the WA is indeed unnecessary, because in the massless 
limit O(a) bulk lattice artifacts are actually absent. 

As for the O(a) boundary terms, we will see that there are two types of 
them: those coming from the action [23] and those coming from the need 
of improving operators inserted at the boundaries, if there is any such local 
factor in O. Boundary terms of the first type will appear in the action with 
coefficients that are even in the Wilson parameter, r. Their non-perturbative 
value is not known. In the standard formulation of the Schrodinger functional 
given by the ^4/paa-Collaboration they have been computed in perturbation 
theory up to two-loops [2D I2H] ■ It turns out that this amount of knowledge 
is numerically adequate for applications. Boundary terms of the second type 
appear multiplied by coefficients that do not have definite r-parity, owing 
to the breaking of parity (see eq. (|3.3|) induced by our choice of fermionic 
boundary conditions (see below). 

As we said, in this paper we will limit ourselves to discuss the mass- 
less theory, which is enough for the non-perturbative evaluation of fermionic 
renormalization constants and the computation of the running of the gauge 
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coupling and quark masses. 

2 The construction of the Schrodinger functional 

The problem with exporting the philosophy of the approach of ref. to the 
standard Schrodinger functional formulation [21] |22j is related to the fact 
that, when the Schrodinger functional is constructed via the iteration of the 
transfer matrix operator |2Z| 1 , the projectors that define which fermionic 
components have to be prescribed at the boundaries are completely deter- 
mined by the form of the transfer matrix and ultimately of the lattice action. 
For instance, for Wilson fermions the boundary conditions must have the 
Dirichlet form 

P + V(x,0) = p(x) , P_V(x,T) = p'(x) , (2.1) 
^(x, 0)P_ = p(x) , ^(x, T)P + = p'(x) , 

with 

p± = i(i±nb). (2.2) 

The choice \r\ = 1 of the Wilson parameter is compulsory, in order for spin 
operators, P±, to be true projectors (Pj = P±). 

Despite the fact that the finite time lattice action defined in this way 
can be shown to enjoy the (spurionic) symmetries that in the infinite volume 
theory were sufficient to guarantee 0(a) improvement of Wilson averages pQ, 
a similar conclusion cannot be drawn here, because upon inverting the sign 
of r the expression of the Schrodinger functional will change by 0(1) terms. 
This operation, in fact, also affects the form of the boundary conditions since 
P-t — > P T under r — > —r. 

The only envisageable way out of this difficulty is to have a formula- 
tion where the structure of the Wilson term and the form of the fermionic 
boundary conditions are not correlated 2 . We then propose to change the 
spin projectors from ()2.2|) and consider a situation where homogeneous TZ^- 
invariant constraints are taken. More precisely, we propose to construct a 
Schrodinger-like functional where fermions are introduced in pairs and obey 
the following homogeneous boundary conditions 

n+^(x,o) = o, n_^(x,T) = 0, 

"^ee ref. I28j for other attempts to define the Schrodinger functional in the continuum 
and on the lattice. 

2 Though developed with the purpose of constructing a lattice Schrodinger functional for 
overlap fermions, the interesting orbifold construction of ref. |29| was of some inspiration 
to us. Ideas similar to the ones discussed here were also presented in ref. 
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V>(x, o)n_ = o , T)n + = o 



(2.3) 



with 



n ± = -(l±r 3 7 5 ) 



(2.4) 



Although not immediately required for 0(a) improvement, we are imagining 
that quarks are introduced in pairs with opposite chiral boundary projectors 
for the two flavour components (say, up and down). As we shall see, this is 
necessary to have a real and positive determinant. 

With only the modification ()2.3)1 of the boundary conditions, the im- 
provement of the expectation value of m.r. 7^-5-even operators (up to 0(a) 
contaminations from boundary lattice artifacts) follows from symmetry ar- 
guments very similar to those that have been used in the infinite volume 
theory. 

We now wish to give a few more details about our construction. Let us 
start by writing down the general form of the action of a pair of massless 
Wilson fermions obeying the boundary conditions ()2.3j) . extended over the 
finite time interval [0, T\. Adapting the analysis of ref. [22] to the present 
situation and separating out the boundary terms from the rest, we write 



Sp — •S'buik + S l B + , 



(2.5) 



T-2a 



S'buik = -yX! H [i>( x ,t)U {x.,t)(r - 7o)V>(x,t + a) + 



x t=a 



+^(x, t + a)(r + 7o)£/ t (x, t)V(*, *) 



+ 



|EE [<M x > *)ffc(x, t)(r - 7 fcMx + k,t) + 



T-a 



x,fc t=a 



+^(x + k, t){r + 7*)tfj£(x, t)*l>(x, t) 



+ 



T-a 



4r 



+a 4 J2 E ^(x, t) M CT (r) + - V(x, t) 

CL 



x £=a 



(2.6) 



S B = -yE [6M x , 0)n + [/ (x, 0)(r - 7o )V(x, a) 



+ 



+^(x, a)(r + 7o)^(x, 0)IT_^(x, 0) 



+ 



x,fc 



^(x)V(x + fc, 0) - K f (x)^(x - fc, 0) , (2.7) 
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S b = -7T E [V>( x > T - a)E7 (x, T — a)(r — 7 o)n+^(x, T) + 
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+^(x, T)n_(r + 7o)E/j(x, T - a)^(x, T 



+ 



+y £ ^(x, T)n_ 7fc [^/(x)^(x + fe,T)- V/ f (x)^(x - k, T)J . (2. 
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x,fc 



For clarity in eqs. (|2.6|) to ()2.8|) we have kept separated spatial vs time com- 
ponents and sums. As usual we are assuming periodicity in space and fix the 
spatial components of the gauge links at the two time boundaries through 



£4(x,0) = ^(x), t4(x,T) = V/(x), fc = l,2,3. (2.9) 



In eq. (12. 6 j) M cr (r) is the critical fermion mass. We remark that no mass 
terms for the fermion fields living at the boundaries have been included. We 
will clarify the reason for that in the next sections. 

With the action ()2.5|) we define a Schrodinger functional through the 
formula 

£W[y/.yi] = jv^ G [U]Jv^^]e- s ^- s \ (2.10) 

where Sym is the pure gauge finite time action j2Ij The gauge integration 
measure, Ufj,c[U], is as explained in refs. |2HE2]- The fermionic integra- 
tion measure, Pjt/p[V>, V>]> is spatially periodic and it is extended in time 
to all fermionic variables from t = a to t = T — a and over the variables 
n_^(x, 0), ^(x, 0)I1 + and II + ^(x, T), ^(x, T)I1_ at the initial and final time, 
respectively. 

It should be noted that this kernel enjoys the usual convolution properties 
only in the continuum limit, unlike the situation one has in the construction 
considered in refs. [2*Tl 12*2*] . Notice that, since the kernel ()2.10|) is not the 
iteration of the transfer matrix, one could even have dropped the terms 
proportional to r in eqs. ()2.7|) and (|2.8|) . 



3 Symmetry properties 

The key observation of this paper is that the action (j2.5j) is invariant under 
a large set of transformations, which leave unaltered the structure of the 
homogeneous fermionic boundary constraints ()2.3|) . They are 

• IZ5 x (r — > — r), where 
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(3.2) 



V : { 



(3.3) 



TZ 5 x V d x V x (V/ <-> V£), where (x D = (-x, T - x )) 

ip(x) — ► e 3i7r//2 ^(a;£i) 
^(x) — > e 3ln / 2 ip(xD) 
U (x) -> Ul{x D - o6) 
t £4(x) C4(z fl - afc) , A; = 1,2, 3 

and P is the standard parity operation (ip = (— x, x )) 

^(x) -> V(^p)7o 
C/ (x) -> C/ (x P ) , 

?7 fc (x) -> C/](x P - aJfe) , A; = 1,2,3 
the product CP, where C is charge conjugation ( T means transposition) 

ip(x) -> i7o72^(a;) T 

V»(a:) -> -^(x) T «7o72 (3.4) 
^(z) E£(a;) , /x = 0,l,2,3 

time inversion around the mid-point T/2, T x (V/ 
(x,T - ar )) 

■0(x) 
■0(x) 

f/o(x) 
U k (x) 



C : 



V£), where (xt 



T 



^(^t)757o 

* Ut(x T - a6) , 

* U k (x T ) , = 



(3.5) 



1,2,3 



• cubic H(3) group 

We explicitely notice that the breaking of V, C and 1Z 5 x V d is entirely due 
to the boundary action terms (|2.7|) and (|2.8|1 . 

Besides the transformations collected above, which are all flavour diago- 
nal, we have also invar iance under 

• the vector rotation in the iso-spin direction 3 

ip(x) -> e iuJT3/2 ijj(x) 
4>(x) -> ^(^)e~^ T3/2 



1a 



(3.6) 



and the four transformations separately acting on only the fermionic 
fields at the boundaries 

(3.7) 
(3.8) 
(3.9) 
(3.10) 



A 3 




• -7 5 r 3 ^(x,0) . 


A 3 


: ^(x,0)- 


-> ^(x,0)7 5 r 3 , 


Al 


: ^(x,T)- 


-> 75 r 3 ^(x,T) , 


A 3 


:^(x,T)- 


> -^(x,T)7 5 r 3 



6 



Finally we mention the reflection symmetry, x (V fc <-> V£), related to 
reflection positivity (see e.g. [H21 IH El)- @ is defined to act as follows 

e{f(u)^(x 1 )...4>(x n )] = 

= r(Q[U})Q[ij(x n )]...e[ij(x 1 )}, (3.11) 
where f(U) is a functional of link variables and for a < xq < T — a 

{9 s /i[ip(x)) = ^(Os/ex)^ 
9 s /e[i>(x)] = joipiOs/ex) 
&s/t[U k (x)] = U£(9 s/e x) 
& s/ AUo(x)} = U^(6 s/e x-a6) 

with 

9 s /t{x, X ) = (x, T - Xq) (3.12) 

to be interpreted as a site (s) or link (£) time-reflection depending on whether 
Nt (T = Ntci) is odd or even, respectively. The action of 6 on the fields at 
the boundaries is 

r e s/ ji_v(x,o)] =^(x,T)n_ 7o 

e s/ ^(x,o)]n + = 7o n + ^(x,r) 

e s/ ,n + ^(x,T)] =^(x,o)n +7o 
{ e s/ ^(x,T)]n_ = 7o n_^(x,o) 

4 Lattice artifacts 

Our aim in this section is to analyze the lattice artifacts affecting the (on- 
shell) expectation value 

(0) T = Jv^Jv^^p^^^e- 3 ™- 3 * , (4.1) 

where O is a m.r. (multi-local) operator. 

The analysis will be based on the assumption that the approach to the 
continuum limit of the Schrodinger functional defined in eq. ()2.10|) can be 
described by a local effective theory. The latter will be characterized a lo- 
cal effective Lagrangian (LEL) including bulk and boundary terms [2Sj- In 
determining the operators contributing to the Symanzik description of the 
lattice expectation value (|4.1|) a key role will be naturally played by the sym- 
metries collected in sect. El Let us now go through the list of the allowed 
operators in the order of increasing dimension. 

• The basic observation is that no dimension 3 operators can be generated 
through radiative corrections in the boundary LEL, because operators of the 
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form ipT (1. / Tb)ifj (with T any of the 16 independent Dirac matrix and with 
or without the insertion of an iso-spin matrix) are all forbidden by some of 
the symmetries above. 

• At dimension 4 there is a number of contributions to the boundary 
LEL. Besides the terms corresponding to the continuum (infinite volume) 
QCD action, there will be 0(a) boundary terms of the form (see eqs. ()2.7|) 
and (EH)) 



The temporal analogs of these terms can be ignored in the present analysis, 
as they can always be traded for the previous ones using the field equations 
of motion. Because of the symmetry TZ§ x (r — > — r) the operators (J4.2)) 
being even under 71$ will intervene multiplied by coefficients that will be 
even functions of r. 

There will also be the 0(a) pure gauge boundary terms 



always with coefficients even in r, while the FF term is excluded by CP. 

Notice that, if in the multi-local operator, O, there are boundary local 
factors, i.e. factors, Ob, which in the continuum limit are kept at vanishing 
distance from the boundaries, there may be (if not excluded by the symme- 
tries listed in sect. EJ) 0(a) terms coming from the collision of Ob with the 
spatially integrated boundary terms of eqs. ()4.2|) and ()4.3|) above. 

• At dimension 5 there are 0(a) effects in the expectation value {0)t 
that in the Symanzik language are described by the insertion, together with 
O, of the integrated (dimension 5 and parity even) bulk LEL density £5. 
Besides these terms there will also appear bulk contact terms, coming from 
the collision under integration of £5 with the various local bulk factors of 
O. In opposition with the previously mentioned local boundary factors, local 
bulk factors are (products of local) operators that in the continuum limit are 
kept at finite distance in physical units from the boundaries. 

Because of the symmetries TI5 x (r — > — r) and 71$ x T>d x V x T, all 
such 0(a) bulk contributions will appear with coefficients odd in r. This 
conclusion follows also here thanks to the fact that in the continuum LEL 
parity is only broken by boundary terms which are not relevant for bulk 
operators, while the bulk LEL density £5 is still even under parity and T and 
consequently odd under 7^5. Thus just like in the infinite volume formulation, 




(4.2) 




(4.3) 
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all 0(a) bulk contributions in the Symanzik expansion will cancel out if the 
expectation value of O computed with r is averaged with its expectation 
value computed with — r, i.e. in the quantity 



(Oh 



WA 



(0) T (r) + (O)r(-r) 



(4.4) 



Notice that in order to compute (0)t|wa it is not necessary to perform two 
independent simulations with opposite value of r. It is enough to notice the 
relation 

(0) T (-r) = (-l) v sM(0) T (r), (4.5) 

where Ps[0] is the parity of O under the transformation TZ 5 . Introducing 
eq. ()4.5|) in (|4.4p . we conclude that for 7^ 5 -even operators 0(a) bulk terms 
are actually absent in the chiral limit and this is why only one simulation is 
needed to get 0(a) bulk improvement. If instead O is an 7?.5-odd operator, 
an identically vanishing result is obtained from the WA. In the absence of 
spontaneous symmetry breaking phenomena affecting the expectation value 
of O, this is indeed the result we would like to get in the continuum. In fact 
(0)t{t) is itself an 0(a) quantity. 



4.1 Boundary operators and O(a) artifacts: an example 

As an example, we briefly discuss a possible setting for the non-perturbative 
computation of the renormalization constant, Zp, of the pseudoscalar density 
operator, ^75 7V0, where b = 1, 2, 3 is an isospin index. The knowledge of Z P 
is relevant, for instance, for the study of the running of the quark mass [2*3*] . 

To extract Zp from simulation data it is sufficient to compute the fol- 
lowing lattice expectation values (for defmiteness we set the isospin index b 
equal to 1) 

($75nVO(0>r/2)$j, 1S2 ) T and ( $ P152 $ P152 >r , (4.6) 

where $p!5 2 and 5>pxS2 are zero three-momentum operators sitting at the 
initial (x = 0) and final (x = T) times, respectively, given by 

$pi52 = « 3 E o)n4( 75 ri + *r 2 )n_^(x, 0) , 

$£152 = a 3 E^(x,T)n4(7 5 ri -2r 2 )n + ^(x,T) . (4.7) 

X ^ 

Notice that $p!s 2 and $p!5 2 are defined in terms of the dynamical boundary 
quark components only. The presence of the projectors H± implies that the 
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Dirac-flavour structure necessarily appears in the combinations (75T1 + ir 2 ) 
and (75T1 — ir 2 ) at x Q = and x = T, respectively. 

An analysis of the dimension three and four boundary operators compati- 
ble with the symmetries listed in sect.EJshows that the fields ()4.7j) renormalize 
multiplicatively and, once summed over their spatial argument, are free from 
0(a) corrections. 

Putting this result together with that on the 0(a) bulk improvement 
reached in sect. EJ one arrives at the conclusion that the calculation of Zp 
carried out with the Schrodiger functional formalism developed in this paper 
is free from 0(a) discretization errors, except of course for those coming from 
the boundary gauge link operators in the first column of eq. ()4.3|) . which 
survive even at vanishing values of the bounday gauge links. 

Although the Schrodinger functional setup proposed in this paper allows 
more freedom (than apparent from the above example) in building fermionic 
boundary operators 3 , the important point we would like to make here is 
that in several cases it may be possible to choose the boundary operators in 
O (those denoted above as Op) in such a way that (0)t is free from 0(a) 
cutoff effects, once the lattice action has been supplemented with boundary 
counterterms. 

5 Fermion determinant 

In this section we want to show that, associated with the fermionic inte- 
gration defined in eq. (|2.1()|) . there is a real and positive determinant. This 
feature is necessary for the proper interpretation of the fermionic contribu- 
tion to the functional integral as a well defined weight for the successive 
gauge integration and the viability of simulations. It is precisely the need 
to fulfill this requirement that has led us to introduce a flavour doublet of 
fermions endowed with opposite chirality boundary conditions 

The situation for the fermion integration is very similar to the one dis- 
cussed in sect. 6 of ref. [22] • Calling u and d the two flavour components 
of if), we can define the two pre-Hilbert spaces Ti u and Hd as the spaces of 
spinors satisfying the conditions dictated by the eqs. (|2.3J) . namely 



With this splitting the action ()2.5p can be rewritten in terms of the operator 



U u = {u I (1 + 75 )u(x, 0) = 0,(1- Ts )«(x, T) = 0} 
H d = {d\ (1- 75 )d(x, 0) = 0, (1 + 75MX, T) = 0} . 



(5.1) 
(5.2) 




(5.3) 



3 This interesting matter is left for a future study. 
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through the formula 4 



S F = (u d)(^ D Q U j ( d u j = uD u u + d~D d d. (5.4) 

One notices that the operator T> admits a well defined eigenvalue problem in 
the Hilbert space 'H u @'Hdi with the functional integration over the fermionic 
variables u and d giving as a result precisely the determinant of D. 

Furthermore we observe that for the purpose of computing the fermionic 
functional integral in (j2.1()j) we can replace the operator (|5.H|) with 

* = ( JU D o" ) ■ < 5 - 5 > 

This replacement only amounts to the harmless change of integration vari- 
ables 

d — > 75G? d — > J75 . (5.6) 

The reason for doing so is that in this way one gets an operator, V, which is 
self-adjoint since 

l5 D dl5 = Dt, (5.7) 

as one can explicitly check with some algebra. Consequently the determinant 
of T> will be a real number. Actually this number, hence the determinant of 
V, is also positive. 

To prove this statement it is convenient to consider the auxiliary operator 

M - ( ° r ° )( ° Du ) - ( ToDi ° ^ (5 

M -\r J I Dl " t d u 1 ' (5 ' 



where T is a block diagonal matrix in the time direction with the following 
structure. All along the diagonal there are Nj- — 1 copies (recall T = NtCl) of 
the 4x4 usual 70 Dirac-matrix, except in the upper-left and the lower-right 
corners where there is a 2 x 2 unit matrix. One checks that Ai maps Ti into 
itself, because 

T D U : Tin — > "H u , 

T Dt : H d ^n d . (5.9) 

Positivity of det[£)] then follows from the chain of equalities 

det[D] = detpD] = det[M] = det[^r ] det[T D u ] = |det[r Dj| 2 , (5.10) 

in which we have used the observation that the determinant of the off- 
diagonal matrix where T intervene in eq. ()5.8|) is equal to 1. 

4 Not to obscure the argument we have suppressed space-time indices. 
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6 Conclusions 



We have shown in this paper that it is possible to modify the usual definition 
of the Schrodinger functional in a way that allows to compute expectation 
values of 7£ 5 -even (multi-local) operators having 0(a) cutoff effects only com- 
ing from the boundaries. Assuming that the latter are under control, either 
because they are just absent, owing to symmetry reasons, or because known 
to some order in perturbation theory, this approach provides a viable scheme 
for the improved computation of the renormalization constants of fermionic 
operators and the running of QCD parameters. 
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